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ABSTRACT 

We discuss evolution of the magnetic flux density and angular velocity in a molecular 
cloud core, on the basis of three-dimensional numerical simulations, in which a rotating 
magnetized cloud fragments and collapses to form a very dense optically thick core of 
> 5 x 10 10 cm" 3 . As the density increases towards the formation of the optically thick core, 
the magnetic flux density and angular velocity converge towards a single relationship between 
the two quantities. If the core is magnetically dominated its magnetic flux density approaches 
1.5(n/5 x 10 10 cm -3 ) 1 / 2 mG, while if the core is rotationally dominated the angular veloc- 
ity approaches 2.57 x 10" 3 (n/5 x 10 10 cm -3 ) 1 ' 2 yr" 1 , where n is the density of the gas. 
We also find that the ratio of the angular velocity to the magnetic flux density remains nearly 
constant until the density exceeds 5 x 10 10 cm" 3 . Fragmentation of the very dense core and 
emergence of outflows from fragments are shown in the subsequent paper. 
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1 INTRODUCTION 

It has long been recognized that magnetic field and rotation affect 
collapse of a molecular cloud, and accordingly, star formation. The 
magnetic and centrifugal forces, as well as the pressure force, op- 
pose the self-gravity of the cloud and delay star formation. Mag- 
netic field and rotation are coupled. Magnetic field is twisted and 
amplified by rotation. The twisted magnetic field brakes cloud ro- 
tation and launches outflows. 

In spite of its importance, only a limited number of numeri- 
cal simulations have been performed for the coupling of magnetic 
field and rotation in a collapsing molecular cloud. The first nu- 
merical simulation o f self- gravitating rotating magnetized clouds 
were performed by iDorfil ll982l) . He found formation of bar- 
like structure for a cloud rotating perpendicular to the magnetic 
field and that of ring-like structure for an aligned rotator. How- 
ever, the grid resolution was limited so that the simulation was 
stopped when the density increased by 200 times from the ini- 
tial value. The spa tial resolution was limited a lso i n other simu - 
lations in 1980's bv lPhillips & Monaghaijjl985t) and lDorfilil989l) . 
who studied the cloud with toroidal magnetic field and that with 
oblique magnetic field, respectively. The spatial resolution was im- 

* E-mail:machida@cfs. chiba-u.ac.jp 
f E-mail:matsu@i. hosei. ac.jp 
| E-mail:tomisaka@th. nao.ac.jp 
§ E-mail:hanawa@cfs. chiba-u.ac.jp 



proved greatly by iTomisakl] 1 19981 l2002h . He considered an ini- 
tially filamentary cloud of n max = 10 4 cm" 3 and followed the 
evolution up to the emergence of magnetically driven outflows 
from the first core of n max > 10 11 cm -3 , where n max denotes 
the maximum density. However, his computation was two dimen- 
sional and could not take account of asymmetry around the rotation 
axis. iBas u & Mouschoviasl J 19941 I l995albl) and iNakamura & Lil 
feOOl 120031) : iLi & Nakamur al <2002l) have got rid of the sym- 
metry around the axis but introduced the thin disk approxima- 
tion. The magnetic braking could not be taken into account in 
these simulations because of the thin disk approximation. Al- 
though Boss <2002l) has performed three-dimensional simulations, 
he has employed approximate magnetohydrodynamical equations. 
The approximation neglects torsion of the magnetic field and re- 
places magnetic tension with the dilution of the gravity. A fully 
three-d imensional nu merical simu l ation h a s just been initiated by 
iMachida. T omisaka & Matsumot cl i2004al) . iHosking & WhitwortM 
(2004), and Matsumot o &Tomisakd<2004l) . 

The recent fully three-dimensional simulations have demon- 
strated that fragmentation of the cloud depends on the magnetic 
field strength. When the magnetic field is weak, a rotating cloud 
fragments after the central density exceeds the critical densi ty, 5 x 
10 10 cm" 3 , i.e., after the formation of Larson's first core iLarsonl 
1969). The magnetic field changes its direction and strength during 
the collapse of the cloud. Thus it is important to study how strong 
a magnetic field the first core has. 
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In this and subsequent papers, we show 144 models 
in which a filamentary cloud collapses to form a magne- 
tized rotating first core. All the models are constructed us- 
in g the fully thre e-dimensi onal numerical s imulation code used 
in iMachida. Tomisaka & Matsumotol <2004al hereafter MTM04). 
This paper shows the evolution by the first core formation stage, 
i.e., the stages before the maximum density reaches the critical 
density, 5 x 10 10 cm -3 . The later stages, i.e., fragmentation of 
the first core and emergence of outfl ows, are shown in the sub- 
sequent paper iMachida. Matsumoto. Hanawa & Tomis aka 2004c, 
hereafter Paperll). 

From analysis of 144 models, we find two variables which 
characterize the evolution of magnetic field and rotation. The first 
one is the ratio of the angular velocity to the magnetic field. This 
remains nearly constant while the maximum density increases from 
5 x 10 2 cm -3 to 5 x 10 10 cm . The second characteristic vari- 
able is the sum of the ratio of the magnetic pressure to the gas pres- 
sure and the square of the angular velocity in units of the freefall 
timescale. This variable converges to a certain value. We refer to the 
convergence as the magnetic flux - spin (B — f2) relation in the fol- 
lowing. We discuss the evolution of the magnetic flux density and 
angular velocity by means of these two characteristic variables. 

This paper is organized as follows: Section 2 denotes the 
framework of our models and the assumptions employed. Section 
3 describes methods of numerical simulations. Section 4 presents 
typical models in the first four subsections and compares various 
models in the last subsection. Section 5 discusses implications of 
the magnetic flux - spin relation and some applications of our mod- 
els to observations. 



2 MODEL 

We consider formation of protostars through fragmentation of a fil- 
amentary molecular cloud by taking account of its magnetic field 
and self-gravity. The magnetic field is assumed to be coupled with 
the gas for simplicity although the molecular gas is only partially 
ionized. Then the dynamics of the cloud are described by the ideal 
magnetohydrodynamical (MHD) equations, 
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where p, v, P, B, and <f> denote the density, velocity, pres- 
sure, magnetic flux density and gravitational potential, respec- 
tively. The ideal MHD approximation is fairly g ood as long as 
the gas density is lower than ~ 10 cm " 3 iNakand [1981 
iNakano. Nishi & Umebava shi 2002). The gas pressure is assumed 
to be 
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where n denotes the number density and is related to the mass den- 
sity p by 



p = 2.3 x 1.67 x 10"" 



(6) 



The critical number density is set to be n cr i = 5 x 10 10 cm 3 
iMasunaea & Inutsukal200fj|) and the sound speed is assumed to be 



c s — 0.19 km s _1 . Thus, this equation of state means that the gas 
is isothermal at T = 10K for n <C n cr i and adiabatic for n 2> n CI i. 

Our initial model is the same as that of Tomisaka (2002) ex- 
cept for the azimuthal perturbation. It is expressed as 
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in the cylindrical coordinates, (r, cp, z) . This initia l model denotes 
a magnetohydrodynamical equilibrium l Stodolkiewicz 1963) when 
Spz(z), 5p v {r, tp) and 8B z (r, <p) are not taken into account. The 
initial density is n Cj o = 5 x 10 2 cm~ 3 on the axis (r — 0). The 
filamentary cloud is supported in part by the magnetic field and 
rotation. This equilibrium is unstable against fragmentation in the 
z-direction. The perturbation in the ^-direction is assumed to be 
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The symbol, A ma x, denotes the wave length of the fastest growing 
perturbation I Matsumo to et al. 1 1 994t) . 

The azimuthal perturbation is assumed to be 



Sp v ,5B,f, — 



A v (r/H) m cos(mv?), for r c < H, 
for r c > H, 



A v cos(mip), 



(14) 



where the azimuthal wavenumber is assumed to be m = 1. The ra- 
dial dependence is chosen so that the density perturbation remains 
regular at the origin (r = 0) at one time step after the initial stage. 
The ratio of density to the magnetic flux density is constant in the 
(^-direction for a given r and z [see equations Q an d OD L 

The initial model is characterized by four nondimensional pa- 
rameters: twice the magnetic-to-thermal pressure ratio, 



a = -B ZCi0 /(47rp c , c Si o), 

the angular velocity normalized by the free-fall timescale, 



(15) 



ui = r2 c , /\/47rGpc,o, (16) 

the amplitude of the perturbation in the ^-direction, A z , and that 
of the non-axisymmetric perturbation, A v . The former two specify 
the equilibrium model, while the later two do the perturbations. We 
made 144 models by combining values listed in TableQ The results 
depend little on the values of A z , thus A z is fixed to be 0.1 in most 
models. 



3 NUMERICAL METHOD 

We employed the same 3D MHD nested grid code as that 
used in MTM04. It incorpor ates the 3D nested grid code of 
iMatsumoto & Hanawal J2003bl) for a hydrodynamical simulation 
and the approximate Riemann solver for the MHD equation 



Collapse and Fragmentation of Molecular Clouds I 3 



jFukuda & Hanawdll999T) . This 3D MHD nested grid code inte- 
grates equations (1) through (5) numerically, with a finite differ- 
ence scheme on the Cartesian coordinates. The solution is a sec- 
ond order accurate, both in space and in time by virtue of the 
Monotone Upstream Scheme for Conservation Law (see e.g., Hirsh 
1990). The Poisson equation is solved by the multigrid iteration 
jMatsumoto & Hanawjl2003al) . We have used Fujitsu VPP 5000, 
vector-parallel supercomputers, for 40 hours to make a typical 
model shown in this paper. 

The nested grid consists of concentric hierarchical rectangu- 
lar subgrids to gain high spatial resolution near the origin. Each 
rectangular grid has the same cell number (= 128 x 128 x 32) 
but a different cell width, h(£) = 2" £ " 5 A max , where I denotes 
the level of the grid and ranges from 1 to ^ max . Thus the coars- 
est rectangular grid of £ = 1 covers the whole computation re- 
gion of — A max < x < A max , —A max < y < A max , and 
z ^ A m ax/2. The solution in z < is constructed from 
that in 2 ^ by the mirror symmetry with respect to z = 0. 
The maximum level number is set at £ m&K — 3 at the initial stage 
(t — 0). A new finer subgrid is generated whenever the minimum 
local Jeans length Aj becomes smaller than /i(^ max )/8. Since the 
density is highest always in the finest subgrid, the generation of the 
new subgrid ens ures the Jeans condition with a margin of a factor 
of 2 (c.f. iTrueloveet al~lll997l) . We have adopted the hyperbolic 
divergence cleaning method of iDedner et al. I {2002) to obtain the 
magnetic field of V ■ B free. 



4 RESULTS 

We have followed all the models shown in this paper until the 
central density exceeds n c > 10 15 cm -3 . This paper describes 
the first half of the evolution for each model, i.e., the stages of 
n c ^ n cr — 5 x 10 10 cm~ 3 . The second half is described in 
the subsequent paper (Paper II). 

Our models are characterized mainly by the strength of the 
magnetic field (a) and the angular velocity (u>). They are classified 
into four groups: (A) models having small a (< 0.1) and small 
ui (< 0.1), (B) those having large a(> 0.1) and small w (< 0.1), 
(C) those having small a (< 0.1) and large u 0.1), and (D) 
those having large a 0.1) and large ui 0.1). In other words, 
the model cloud has a weak magnetic field and rotates slowly in 
group A, while it has a relatively strong magnetic field and rela- 
tively large angular momentum in group D. Each group is described 
separately in the following subsections, in each of which two typi- 
cal models of A v = 0.01 (S) and 0.2 (L) are shown. The typical 
models are named after the group (A, B, C, or D) and A v (S or L). 
Model AS has a = 0.01, u) = 0.01, and A v = 0.01, for example. 
Table |2| shows the values of a, lo, A z , and A v for the 8 typical 
models shown in the following subsections. It also shows the initial 
magnetic field (B ZCj o), the initial angular velocity (fi c o), the wave- 
length of the perturbation in the ^-direction (A max ), the mass (M) 
of the gas contained in the region of \z\ ^ A max /2, and the epoch 
at which the density becomes infinity (t{). 



4.1 Weak Magnetized and Slowly Rotating Cloud 

This subsection displays model AS as a typical model having a 
weak magnetic field and slow rotation. Model AS has parameters 
a = 0.01, uj = 0.01 and A v = 0.01. Fig. □ shows the cloud 
evolution in model AS by a series of cross sections. 



As shown in Fig. Q a gas cloud is transformed from a pro- 
late one to an oblate one on the y = plane (see lower panels), 
while it maintains a round shape on the z — plane (upper pan- 
els) in the period of 5.5 x 10 2 cm -3 < n c < 7.6 x 10 10 cm -3 . 
The velocity field is almost spherically symmetric while the central 
density increases from 5 x 10 5 cm" 3 to 2 x 10 9 cm . The dense 
cloud is prolate and elongated in the z-direction in the lower panel 
of Fig. Q(b), while it is nearly spherical in the lower panel of Fig.Q 
(c). In this early collapse phase, the cloud contracts along the ma- 
jor axis (i.e. z-axis), regardless of the magnetic field and rotation as 
discussed in lBonnell et al. I ll99cl) . An oblate core is seen in panel 
Fig. Q(e) and a thin disk is seen in the lower panel of Fig. Q(f). 
The collapse is dynamical at the stages shown in Fig.Q(c) through 
Fig. 0(f). The radial infall velocity reaches v r = —0.52 km s~ 
on the 2 = plane in Fig.Q(e), while the vertical infall velocity 
does v z = ±0.58 km s" 1 on the z-axis. The rotation velocity is 
v v = 0.047 km s" 1 at maximum and much smaller than the infall 
velocities. This means gas contracts spherically in this phase. The 
difference between the radial and vertical infall velocities is still 
small (|« r ,max| = 0.61 km s" 1 and ju z , max | = 0.8 km s" 1 ) in 
Fig.Q(f) although a high density disk is formed. 

The density increase is well approximated by p c = 
2.2/[47rG(t - if) 2 ] in the period of 5 x 10 4 cm" 3 sC n < 
10 9 cm" 3 as shown by the thick solid curve in Fig.|2l The offset 
is taken to be tf = 5.96 x 10 6 yr so that the central density in- 
creases in proportion to the inverse squar e of the time in the widest 
span in log p c , as shown in Larson 1 1969). Remember that the sim- 
ilarity solution of Larson (1969) and Penston (1969) gives p c — 
1.667/[4-7rG (t — tf) 2 ] for spherical collapse of anon-magnetized 
non-rotating isothermal cloud. We have checked that the density 
increase is well approximated by p c ~ 1.65/[47rG(t — tf) 2 ] in a 
non-magnetized and non-rotating cloud of our test calculation. The 
density increase is 15 % slower in model AS than in the similarity 
solution, since (2.2/1.667) 1 / 2 ~ 1.15. This small difference is 
due to the rotation and magnetic field. 

To evaluate the change in the core shape shown in Fig. Q 
we measure the moment of inertia for the high density gas 
of p ^ 0.1p c . We derive the major axis (hi), minor axis 
(h„), and z-axis (h z ) from the moment of inertia according to 
Matsumoto & Hanawa ( 1999). The oblateness is defined as e b = 
{hihs) 1 / 2 /h z and the axis ratio is defined as e ar — hi/h s — 1. 

The oblateness is denoted by the thick solid curve as a function 
of time in Fig.|5|(a) and as a function of the central density in Fig.|5| 
(b). The oblateness is nearly constant at £ oc . = 0.27 in the period 
oft < 4 x 10 6 yr (or 5 x 10 2 cm" 3 < n c < 5 x 10 3 cm" 3 ). It 
increases and reaches e b — 1 at the stage of n c — 2 x 10 6 cm" 3 , 
which is shown in the lower panel of Fig. 0(c). The oblateness 
reaches e Q b = 2.9, when the disk-like structure is formed at n c = 
7.6 x 10 10 cm" 3 as shown in the lower panel of Fig. Q(f). The 
increase in e D b is monotonic in the period of n c > 5 x 10 3 cm" 3 . 

The axis ratio is denoted by the thick solid curve as a function 
of time in Fig. [3] (c) and as a function of the central density in 
Fig.|3|(d). The axis ratio decreases from e ar = 0.01 to 7 x 10" 4 
after oscillating once over the period of t ^ 5 x 10 6 yr ( or n c < 
5 x 10 3 cm" 3 ). It increases in proportion to over the period 
of n c > 5 x 10 3 cm" 3 . The growth rate of e ar coincides with 
that of the bar mode growing in the spherical runaway collapse 
(Hanawa & Matsumoto 1999). The axis ratio grows up to 3.5 x 
10" 3 in the isothermal collapse phase as shown in Fig.|3](d). In 
order to examine dependence on the axis ratio, we compare models 
AS and AL, of which model parameters are the same except for the 
amplitude of the non-axisymmetric perturbation, A v . The value of 
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e ar is 20 times larger in model AL than in model AS at a given 
stage and reaches e ar = 6.8 x 10" 2 at n c — n C ri- The axis ratio 
is proportional to A v . The oblateness is nearly the same in models 
AS and AL. The non-axisymmetric perturbation grows linearly in 
models AS and AL. 

The magnetic flux density increases as the density increases. 
The left panel of Fig. |4| shows the square root of the ratio of the 
magnetic pressure to the gas pressure, B zc / (87rp c c 2 ) 1//2 , as a func- 
tion of n c . Note that the ordinate is normalized by the initial value. 
It increases in proportion to one sixth the power of the density, i.e., 
B zc /(87rpcC 2 ) 1,/2 oc nj^ 6 , in the period of 10 6 cm -3 ^ n c ^ 
10 9 cm~ 3 . This means that the magnetic field increases in pro- 
portion to B zc oc n 2,/3 . This increase in _B ZC is consistent with 
the spherical collapse of the core. When the collapse is spherically 
symmetric, the density and magnetic field increase inversely pro- 
portional to the cubic and square of the radius, respectively, since 
the magnetic field is frozen in the gas. Hence, the magnetic field is 
proportional to two thirds the power of the density, B zc oc fij 1 '. 

After the central density exceeds 10 9 cm~ 3 , the growth of 
the magnetic field slows down. This slowdown coincides with the 
change in the core shape. The core is significantly oblate in the 
period of n c > 10 9 cm~ 3 . Remember that the magnetic field is 

1/2 

proportional to the square root of the density (B zc oc n c ) when 
a magnetized disklike gas cloud collapses (Scott & Black 1980). 
This is because the disk is nearly in a hydrostatic equilibrium in the 
z-direction and the isothermal disk has the relation n oc E 2 . We 
use the terminology, the "disk collapse", for this radial collapse of 
a disklike gas cloud. In the disk collapse, the magnetic flux density 
increases in proportion to the surface density (_B ZC oc E c ) since the 
gas is frozen in a magnetic flux tube. The relations, n c oc E 2 and 
B zc oc E c , yield B zc oc nj' 2 . In the period of n c > 10 9 cm -3 , the 
growth rate of the magnetic field is intermediate between those for 
the spherical collapse and for the disk collapse. This is consistent 
with the density change over the same period. 

As well as the magnetic flux density, the angular velocity 
of the core increases as the density increases. The right panel of 
Fig. 4 shows the ratio of the angular velocity to the magnetic field 
(Q c /B zc ) normalized by the initial value (fl c ,o/B zc ,o) as a func- 
tion of n c . The ratio is nearly constant at the initial value. This is 
because both the specific angular momentum (j) and the magnetic 
flux (<J>) are conserved for a central magnetic flux tube. Both the 
angular velocity and magnetic field increase proportionally to the 
inverse square of the tube radius. Hence the ratio is constant in both 
the spherical and the disk collapse. The conservation of the specific 
angular momentum implies that none of the magnetic torque, grav- 
itational torque, and (/^-component of the pressure force are signifi- 
cant. 

Since fl c /B zc is nearly constant, the angular velocity in- 
creases in proportion to n 2 ^ 3 in the period of 10 6 cm -3 ^ n c ^ 
10 9 cm -3 and the growth of f2 c slows down in the period of 
n c > 10 9 cm" 3 . When measured in units of the free fall timescale, 
the angular velocity increases in proportion to f2 c (47rGpc)~ 1,/2 oc 
nd % in the former period. The angular velocity in units of the free 
fall timescale denotes the square root of the ratio of the centrifu- 
gal force to the gravitational force. The magnetic field and rotation 
strengthen in the same manner during the spherical collapse, since 
both B zc (8np c cl)~ 1 ^ 2 and f2 c (47rGp c ) _1 ' /2 increase in propor- 

1/6 

tion to p c ■ 

Model AS is similar to model B of iMatsumoto et at. I J 19971) . 
although our model AS includes a very weak magnetic field. The 
magnetic field influences little the cloud collapse. Fig.|5|shows that 



the magnetic field is not twisted but pinched at the stages of n c 
10 B cm -3 as shown in panels (d) - (f) of Fig. [5] Each panel denotes 
the magnetic field lines for the corresponding stage shown in each 
panel of Fig.Q This weak magnetic field has no significant effect. 
When a < 0.1 and u < 0.1, the effects of magnetic field and 
rotation are very small. 



4.2 Strongly Magnetized and Slowly Rotating Cloud 

Model BL is shown as a typical example of models in this sub- 
section having large a and small uj. Model BL has parameters 
a — 0.1, oj — 0.01 and A v = 0.2. The parameters of model BL 
are the same as those of model AL except for a, which is 0.01 for 
model AL and 0.1 for BL (Table|2j. When a ^ 0.1, the magnetic 
pressure becomes comparable to the gas pressure in the course of 
cloud collapse and decelerates the radial collapse significantly. The 
magnetic braking is also effective in models BL and BS. 

Also in model BL the high density core changes its form from 
prolate to oblate as the central density increases, as shown in Fig.|6| 
which is the same as Fig.Qbut for model BL. The change in the 
core shape is due to the magnetic field, which is amplified during 
the spherical collapse. The ratio of the magnetic pressure to the gas 
pressure is 0.1 1 at the stage of n c = 2 x 10 6 , while it is only 0.05 
at the initial stage. Each panel of Fig. [6] denotes the density and 
velocity distribution at the stage of (a) n c = 8.2 x 10 3 cm -3 , (b) 
5.6 x 10 4 cm" 3 , (c) 7.9 x 10 6 cm~ 3 , and (d) 6.0 x 10 10 cm -3 . 
At the stage of n c — 8.2 x 10 3 cm~ 3 , the oblateness is e D b = 
0.58 in model BL [Fig. 6(a)] while e ob = 0.45 in model AL [Fig. 
1(b)] . The core is more oblate in model BL than in models AL and 
AS when compared at a given stage with the same central density 
[Fig. 3(b)]. The oblateness is e Q b = 5.3 at the stage of n c = 5 x 
10 10 cm" 3 in model BL, while e ob = 2.9 in model AS (see Fig. 3). 
The oblateness increases slowly over the period of 3 x 10 3 cm -3 < 
n c < 10 10 cm -3 in model BL and is saturated around e b — 5 in 
the period of n c > 10 10 cm -3 . 

As well as in models AS and AL, the axis ratio decreases 
from the initial value of e ar = 0.2 to 0.015 in the period of 
n c 10 4 cm" 3 in model BL. Then it switches to growing in 
proportion to n\^' in the period of n c > 10 4 cm -3 . The core is 
elliptic on the x — y plane and the axis ratio is e ar = 0.23 at the stage 
of n c — 5 x 10 10 cm -3 , as shown in the upper panel of Fig.|6|(d). 
The amplitude of the non-axisymmetic perturbation is linearly pro- 
portional to the initial amplitude. The axis ratio is always smaller 
by a factor of 20 in model BS than in model BL when compared 
at the stage of a given central density. Models BS and BL have the 
same model parameters except for A v . 

Since the core is appreciably oblate, the infall velocity is 
higher in the vertical direction than in the radial direction. At 
the stage of n c = 5 x 10 10 cm" 3 the radial infall velocity is 
v r = —0.46 km s" 1 at maximum in z — plane while the verti- 
cal infall velocity is v z — 0.7km s" 1 at maximum on the z-axis. 
The radial infall velocity is smaller and the vertical infall velocity 
is larger than in model AS. This asymmetry is due to the magnetic 
field. The rotation velocity is quite small (v v = 0.03 km s" 1 ) and 
the centrifugal force is negligible. The density increase due to col- 
lapse is slower in model BL than in model AS. The growth of the 
central density is well approximated by p c = 2.5/[47rG (t — U ) 2 ]. 
The growth rate is 7 % smaller than that of model AS at a given 
central density. 

Also in model BL, the magnetic field strengthens as the den- 
sity increases. The ratio of the magnetic pressure to the gas pressure 
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increases very slowly in the period of n c < 10 6 cm -3 (see Fig. 4). 
It is saturated around (B zc / B ZCt o) 2 /(p c /pc,o) — 1.3 in the pe- 
riod of n c > 10 6 cm -3 . The magnetic field decelerates the radial 
collapse appreciably as shown earlier. 

Fig. shows the magnetic field lines at the stage of n c = 
6 x 10 10 cm -3 . The magnetic field lines break at the levels of 
z ~ 15 AU and 30 AU near the disk surface. The latter break 
corresponds to a fast-mode MH D shock, which is ess entially the 
same as the shock waves seen in Norman. Wilson & Barton 1 1980), 
iMatsumoto et al. I Jl997t) . and iNakamura et al. I Jl999l) . They are 
squeezed and vertical to the midplane below the shock front while 
they are open above the shock front. The disk formation is due to 
the magnetic field. 

The magnetic field extracts angular momentum from the core. 
As shown in Fig.[4](b), the ratio of angular velocity to the magnetic 
field decreases by 30 % in the period of 5 x 10 2 cm -3 < n c 
5 x 10 10 cm" 3 in model BL, although it remains constant in model 
AL. The decrease is due to the magnetic braking. The twisted mag- 
netic field transfers the angular momentum of the core outwards. 
The specific angular momentum of the core is 70 % of the initial 
value at the stage of n c = 5 x 10 10 cm -3 . The angular velocity nor- 
malized by free-fall timescale [Q, c / (47rGp c ) 1//2 ] increases slightly 
from 0.01 to 0.015 in the period of 5 x 10 2 < n c < 5 x 10 10 cm" 3 
in model BL, while it spins up from 0.01 to 0.06 in model AL. We 
discuss this difference again in §4.5 in which we compare the in- 
crease in Q. c for various models. 

The efficiency of the magnetic braking is qualitatively similar 
in models BS and BL. 

Models BL and BS are similar to model C o flTomisakiJiT99l 
and model B 1 of Nakamura et al. 1 <1999h . although the earlier mod- 
els include neither rotation nor non-axisymmetric perturbation. The 
rotation and non-axisymmetric perturbation have little effect on the 
cloud collapse in our models BL and BS. When the initial mag- 
netic pressure is larger than a tenth of the gas pressure (a ^ 0.1), 
initially weak magnetic field is amplified during the collapse and 
affects the evolution of the core. The magnetic pressure decelerates 
the radial collapse and leads to disk formation. Also the magnetic 
braking is appreciable. 

4.3 Weakly Magnetized and Rapidly Rotating Cloud 

This subsection describes model CS as a typical example for mod- 
els having small a and large u). Model CS has parameters of 
a = 0.01, to = 0.5, and A v = 0.01. When uj > 0.1, rota- 
tion affects the collapse of the cloud significantly. 

In model CS, a rotating disk forms at an early stage of low 
central density. Each panel of Fig. [8] denotes the density and ve- 
locity distribution at the stages of (a) n c = 5.2 x 10 3 cm' 3 , (b) 
6.5 x 10 4 cm" 3 , (c) 5.7 x 10 6 cm" 3 , and (d) 8.3 x 10 10 cm -3 . 
The rotating disk is clearly seen at the stage of 6.5 x 10 4 cm -3 . The 
oblateness reaches e D b = 3.0 at the stage of 4.0 x 10 4 cm" 3 and 
is saturated around e D b — 3.5 in the period of 5 x 10 5 cm -3 < 
n c < 5 x 10 10 cm -3 as shown in Fig.[3](b). 

The axis ratio increases up to e ar = 0.2 by the stage of n c = 
5 x 10 10 cm -3 in model CS (see Fig. 3). At the stage of n c = 
5 x 10 10 cm -3 , the axis ratio is larger in model CS than in model 
AS, while it is the same at the initial stage. The difference arises in 
the period of n c ^ 5 x 10 4 cm -3 . The axis ratio remains around 
0.01 in model CS, while it decreases to 7 x 10" 4 in model AS. 
The axis ratio grows roughly in proportion to n~J 6 in the period of 
n > 5 x 10 4 cm" 3 both in model AS and in model BS. We have 
confirmed that the non-axisymmetric perturbation is proportional 



to the initial perturbation by comparing with model CL of which 
initial parameters are the same as those of model CS except for 
A v . The axis ratio is 20 times larger in model CL than in model 
CS in the period of n c < 1.0 x 10 cm -3 . The axis ratio reaches 
Ear = 10.2 and the high density core has a bar shape at the stage of 
h c = 5x 10 10 cm" 3 in model CL. 

The increase in the central density is approximated by p c = 
6.2/[4ttG (t - if) 2 ]. The rate of the increase is appreciably smaller 
than those of models AS and BS. It is 1.93 times smaller than that 
of the spherical collapse at a given central density. The relatively 
slow collapse is due to fast rotation. 

In the period of n c < 5 x 10 5 cm" 3 the cloud collapses mainly 
in the vertical direction along the magnetic field. Accordingly the 
magnetic field increases a little and the ratio of the magnetic pres- 
sure to the gas pressure decreases in this period (see Fig.[4j- Note 
that the square root of the ratio of the magnetic pressure to the 

— 1/2 

gas pressure decreases in proportion to p c when the collapse is 
purely vertical along the magnetic field. Also, the angular velocity 
increases a little and decreases in proportion to pc 1 ^ 2 when mea- 
sured in the freefall timescale. 

In the period of 5 x 10 2 cm" 3 < n c < 5 x 10 10 cm" 3 , 
the magnetic field (B zc ) strengthens and the angular velocity of 
the core (fi c ) continue to increases. However, the ratio of the 
magnetic pressure to the gas pressure remains nearly constant 
around (B zc /B zc ,o) 2 /{fie/ Pc,o) = 0.5 for 5 x 10 5 cm" 3 < 
n c < 5 x 10 10 cm" 3 as shown in Fig. [4] The angular velocity 
measured in the freefall timescale is also nearly constant around 
r2 c (47rGp c )~ 1//2 = 0.2. In other words, both the magnetic field 

1/2 

and angular velocity increase in proportion to p c ■ These depen- 
dences of B zc and Q c on p c indicate that the core collapses in the 
radial direction while maintaining a disk shape. They are the same 
as those in the similarity solution for a self-gravitationally collaps- 
ing gas disk (Tomisaka 1995, 2002; Nakamura et al. 1995; Mat- 
sumoto et al. 1997; Saigo & Hanawa 1998). 

The ratio of the angular velocity to the magnetic field is con- 
stant in the period of n c < 10 7 cm" 3 . It increases up to 1.2 by 
the stage of n c = 5 x 10 10 cm" 3 . This increase is due to the 
torsional Alfven wave. The magnetic braking is not significant in 
model CS. The ratio of the angular velocity to the magnetic field is 
constant during the collapse as shown in Fig.|4](b). This confirms 
that the specific angular momentum is conserved. The magnetic 
field is twisted by fast rotation as shown in Fig.|9| It is also bent at 
the shock front as well as in model BS. The twisted magnetic field 
is too weak to have any appreciable dynamical effects. 

The infall velocity is higher vertically than radially. The max- 
imum infall velocity is v r ,max = —0.37 km s" 1 radially and 
^z.max = ±0.59 km s" 1 at the stage of n c — 5 x 10 10 cm" 3 . 
The maximum rotation velocity is v v = 0.36 km s" 1 and exceeds 
the sound speed at the same stage. Thus both the infall and rotation 
are supersonic. This dynamically infalli ng gas disk is si milar to in- 
falling envelo pes observed in HL Tau ( Havashi et al. 1 119931) and 
L1551 IRS5 i Ohash i et al. Ill996t ISaito et al. 1 11996) in the sense 
that the radial infall velocity is comparable with the rotation veloc- 
ity. The vertical inflow along the z-axis forms shock waves twice, 
once at the stage of n c = 5.7 x 10 6 cm" 3 [see Fig.[8|(c)] and at 
that of n c = 8.3 x 10 10 cm" 3 [see Fig.|S|(d)]. The former forms 
at z = ±4 x 10 3 AU, and the latter at 2 = ± 40 AU. These 
shock waves are essentially the same as those seen in Matsumoto 
et al. (1997). The oblateness has a temporal maximum value at the 
stages of the shock formation. 
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4.4 Strongly Magnetized and Rapidly Rotating Cloud 

This subsection describes model DL as a typical example of mod- 
els having large a and large uj. Model DL has parameters a — 1.0, 
uj = 0.5 and A v = 0.2. When a ^ 0.1 and w ^ 0.1, both mag- 
netic field and rotation affect the collapse of the cloud significantly. 
The magnetic braking is also effective. 

Fig. llOl shows formation of a magnetized rotating disk that de- 
forms to an elongated high density bar in model DL. Each panel 
denotes the density and velocity distribution at the stage of (a) 
n c = 5.6 x 10 3 cm" 3 , (b) 6.8 x 10 4 cm" 3 , (c) 5.1 x 10 6 cm" 3 , 
and (d) 5.3 x 10 10 cm -3 . The high density gas has an oblateness 
of e b = 4.2 at the stage of n c = 6.8 x 10 4 cm -3 . The oblateness 
reaches its maximum at n c ~ 4 x 10° cm" 3 and oscillates around 
e ob ~ 5 in the period of 4 x 10 5 cm" 3 < n c < 5 x 10 10 cm" 3 
(see Fig.|5J. As a result of the strong magnetic field and fast rota- 
tion, the disk forms at an earlier stage in model DL than in the other 
models shown in the previous subsections. 

The increase in the central density is approximated by n c = 
4.9/|4tt G (t - tf) 2 \. This is slower than in models AS and BS, 
although faster than in model CS. The density increase is faster in 
a model having a larger a for a given uj ^ 0.1. This is because 
a stronger magnetic field brakes the rotating core more effectively 
and the centrifugal force is reduced more. Remember that the den- 
sity increase is slower in a model having a larger a, when u < 0.1. 
When the angular momentum of the cloud is very small, the cen- 
trifugal force is negligible and its reduction due to the magnetic 
braking is unimportant. A stronger magnetic field decelerates the 
collapse through higher magnetic pressure and tension. The mag- 
netic field plays two roles: acceleration of the collapse through 
magnetic braking, and deceleration of the collapse through mag- 
netic pressure and tension. The former dominates for uj ^ 0.1 
while the latter dominates for uj < 0.1. 

As shown in Fig. 1101 (d) upper panel, the disk is elongated 
into a bar of e ar = 15 at the stage of 5.3 x 10 10 cm -3 . As well 
as in model CL, the axis ratio remains nearly constant at the be- 
ginning and increases in proportion to n\ from an early stage 
of n c < 5 x 10 3 cm -3 in model DL (see Fig.[3J. An elongated 
bar forms in the models in which the non-axisymmetric perturba- 
tion is relatively large {A v > 0.2) and does not diminish in the 
early phase. Remember that the axis ratio decreases in the period 
of n c < 1-0 x 10 4 cm -3 in models AL and BL. The axis ratio in- 
creases in proportion to nV & in all models while the core collapses 
dynamically. The final axis ratio depends on the initial value and 
the amount of damping in the early phase. The initial damping is 
smaller when either the initial magnetic fiel d or rotation is larger. 

A similar bar structure is also seen in iDurisen et ai~l fl9861 
and Bate 1 1998). The bar structure develops as a result of the bar 
mode instability, when the ratio of rotational to gravitational en- 
ergy (/3) of the core exceeds /3 = 0.274. This /3 is related to uj by 
f3 = uj 2 , when the cloud is spherical and has constant density and 
angular velocity. Thus, the criterion for the 'bar mode instability', 
> 0.274, corresponds to u) > 0.523. The value of uj continues to 
decrease until it converges to uj ~ 0.2 as denoted in the following 
subsection. Therefore, the condition uj > 0.523 is never fulfilled, 
and the bar should form in mode l D by another mechanism . While 
the bar mode instability of lHanawa & Matsumoto! ^999) works 
o nly in a dy namically collapsing cloud , the bar mode instability 
of lDurisen etall 1 1986) and lBatel il993l does in a cloud in hydro- 
static equilibrium. 

There is another evidence that the bar is formed not by fast 
rotation in our models. The bar forms also in the model of (a, w, 



A v ) = (3, 0, 0.2), which is listed as model 56 in Table 2 of Paper II. 
The bar formation can not be due to rotation since the cloud does 
not rotate in this model. (See Table 2 of Paper II for the list of the 
models in which the bar forms at the end of the isothermal phase.) 

As well as in models CS and CL, the vertical infall dominates 
over the radial infall in the period of n c < 7 x 10 4 cm" 3 in models 
DS and DL. The ratio of the magnetic pressure to the gas pres- 
sure normalized by its initial value decreases to 0.5 in the period 
as shown in Fig.|4| Then it oscillates around 0.5 in the period of 
10 5 < n c < 5 x 10 10 cm . The epoch of disk formation coin- 
cides with that at which the ratio of the magnetic pressure to the 
gas pressure reaches its first local minimum value. 

The vertical inflow also forms shock waves twice in model 
DL. Fig.Qotc) shows the outer shock located at z = ±7000 AU. 
The flow is nearly vertical above the front while it is horizontal 
below. The epoch of shock formation coincides with that of the 
temporarily maximum oblateness as in model CL. 

The magnetic braking slows the spin of the collaps- 
ing disk in model DL. The initial central angular velocity is 
n c /(47rGp c ) 1/2 = 0.5 in both models CL and DL. The central 
angular velocity decreases to n c /(47rGp c ) 1 ^ 2 = 0.21 by the stage 
of n c = 5 x 10 4 cm" 3 in model DL, while it decrease to 0.24 in 
model CL. The ratio of the angular velocity to the magnetic field 
(Q, c /B zc ) normalized by the initial value (Q c ,o / B ZCj o) decreases 
to 70 % of the initial value in model DL (see Fig. [4}- The mag- 
netic braking is effective in the period of n c ^ 7 x 10 5 cm -3 as 
in model BL. The ratio of the angular velocity to the magnetic field 
increases in the period of 7 x 10 5 cm" 3 < n c < 5 x 10 8 cm -3 . 
This spin is due to the magnetic torque. The twist of the magnetic 
field is bounded by the shock front and the torsional Alfven wave 
is reflected there. Thus the angular momentum is not released from 
the core in model DL. 

The ratio of the magnetic pressure to gas pressure decreases 
from 0.5 to 0.1 in the period of 5 x 10 2 cm" 3 < n c < 7 x 
10 s cm -3 , and remains around 0.1 in the period of n c > 7 x 
10 s cm -3 . Thus, importance of the magnetic force relative to the 
centrifugal force decreases in models DS and DL. Fig. ll ll illustrates 
the magnetic field for the stages shown in Fig. 1101 The magnetic 
field lines are twisted but less pinched than in model CS. They are 
twisted at a higher z in model DL than in model CS. As shown in 
Fig. ll 11 (d). the magnetic field is squeezed to stem vertically from 
the bar and the magnetic flux density is large in the bar. In Fig. ll II 
(d), magnetic field lines are bent at z ~ 40AU, which corresponds 
to the shock front. Inside the shock front, the magnetic field lines 
are ran vertically and hardly twisted, while twisted moderately out- 
side of the shock front. 

Models DL and DS have the same initial model parameters 
except for A v , which is smaller by a factor of 20 in model DS. 
As a result, only the axis ratio differs appreciably between models 
DL and DS. A high density disk is seen at the stage of n c = 5 x 
10 10 cm -3 in model DS while the elongated bar is seen in model 
DL. The axis ratio is smaller by a factor 20 in model DL throughout 
the evolution. 

4.5 Magnetic Flux - Spin Relation 

The filamentary cloud fragments to form a high density core of 
n c > 5 x 10 10 in all the models computed. The formation of the 
core is dynamical and the central density increases in proportion 
to the inverse square of the time. As the central density increases, 
the magnetic field increases roughly in proportion to a power of p c . 
The power index, however, differs and depends on the geometry of 
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the collapse as shown in the previous subsections. Also the angu- 
lar velocity increases in proportion to a power of p and the power 
index depends on the geometry of the collapse and the strength of 
the magnetic field. To summarize the increase in B zc and £l c we 
have plotted the evolutionary locus of the core in Fig. 12. The ab- 
scissa denotes the square root of the ratio of the magnetic pressure 
to the gas pressure, B zc (87rCap c ) _1//2 , in the logarithmic scale. The 
ordinate denotes the angular velocity normalized by the freefall 
timescale, fi c (47rGp c ) _1 ^ 2 , on the logarithmic scale. Each curve 
denotes the evolutionary locus for a model. The asterisks denote the 
initial stages. The circles, squares and triangles denote the stages of 
n c = 5 x 10 4 cm -3 , 5 x 10 6 cm -3 , and 5 x 10 s cm -3 , respec- 
tively. Models without magnetic field are shown inside the upper 
left box. Models without rotation are shown inside the lower right 
box. 

The evolutionary loci are systematically ordered in Fig. 12. 
They are aligned to evolve toward the upper right in the lower left 
region. Models AL and AS belong to this region of weak magnetic 
field and slow rotation. On the other hand, the loci are aligned to 
evolve toward the lower left in the upper region (fast rotation) and 
in the right region (strong magnetic field). Models CL, CS, DL, 
and DS belong to these regions. Models BL and BS appear in the 
middle of the panel. Their loci are nearly horizontal and the angular 
velocity measured on the freefall timescale does not increase as a 
result of the magnetic braking. 

We can deduce several rules for the collapse of a magnetized 
rotating gas cloud from Fig. 12. First all the loci seem to converge 
on the curve, 



fi 



+ 



B. 



(17) 



(0.2) 2 4nGp c (0.36) 2 8nc 2 8 p c 

Equation J 171 is denoted by the thick solid curve in Fig. 12. We 
call this curve the magnetic flux - spin relation or B — Q. relation in 
the following. The first term of equation Ml\ is proportional to the 
square of the angular velocity normalized by the freefall timescale 
and accordingly is proportional to the ratio of the centrifugal force 
to the gravity. The second term is proportional to the ratio of the 
magnetic pressure to the gas pressure. The numerators are propor- 
tional to the anisotropic forces which suppress only the radial infall. 
Convergence to equation < 171 indicates that the sum of the centrifu- 
gal and magnetic forces are regulated to be at a certain value. This 
rule involves the rules found by Matsumoto et al. (1997) and Naka- 
mura et al. (1999) as a special case. The former showed that the ra- 
tio of the specific angular momentum to the core mass converges to 
a half of the critical value j = 0.5(27rGAf/c s ) for models having 
no magnetic field. The latter showed that the ratio of the magnetic 
field to the surface density tends to be a half of the critical one, 
i.e., -B zc = 0.5(27rG) 1//2 E, for collapse of a non-rotating cloud. 
See the models shown inside the upper left box and those inside the 
lower right box to confirm that they also converge to equation I17i . 

The magnetic flux - spin relation is related to formation of the 
shock waves. The first shock wave forms exactly when the evolu- 
tionary locus reaches the B~Q relation. After the shock formation, 
the evolutionary locus leaves it temporarily and reaches it again at 
the formation of the second shock wave. The shock strength is also 
related to the distance between the initial stage and the B — Q. rela- 
tion on the B — Q. diagram. The shock wave is stronger in a model 
starting from a more distant place from equation 1171 on the dia- 
gram. No shock wave forms in a model of which the initial stage 
is close to equation <17l on the diagram (see, e.g., the model of 
a — 0.1 and u = 0.1 shown in Fig. 12). 

Second, the magnetic braking is appreciable only in mod- 



els having a ~ 0.1. The effect of the magnetic braking is 
evaluated from the slope on the diagram, dlog £l c /d log B zc . 
When the specific angular momentum is conserved, the slope is 
dlog Q c /d\og Bz C — 1 as discussed in subsection 4.3. The slope 
differs appreciably from unity near the lower part of the B — Q 
relation in Fig. 12. It is appreciably smaller than unity on the left- 
hand side of the B — Q relation, while it is appreciably larger on 
the right-hand side. When the initial magnetic field is very weak, 
its magnetic torque is negligible. When the initial magnetic field is 
strong, the vertical collapse dominates. The magnetic braking re- 
duces the specific angular momentum by 30 — 40% by the stage of 
n c = 10 6 cm -3 . However, it does not operate effectively beyond 
the stage. We will discuss the implication of equation I17i in the 
next section. 



5 DISCUSSION 

As shown in the previous section, the magnetic flux density and 
angular velocity converge on equation <17l in Fig. 12, the B — 
fi diagram. Thus, we can evaluate the magnetic flux density and 
angular velocity of the first core to be 



2.57 x 10" 3 yr- 1 



B z , 



1.50 x 10 4 ^G 



= 1, 



(18) 



by substituting p c — 1.92 x 10 13 gem 3 (equivalent to n c = 
5 x 10 10 cm -3 ) and c s = 0.19 km s~\ into equation IT71 . 
equation 1 1 8i implies that the first core has either the 'standard' 
magnetic flux density (15 mG) or the 'standard' angular veloc- 
ity (2.57 x 10 -3 yr _1 ), unless the initial magnetic field is very 
weak and the rotation is very slow. When both the magnetic flux 
density and angular velocity are negligible, the cloud collapses al- 
most spherically and hence both B zc and fi c increase in propor- 

2/3 

tion to n c . Thus, either the magnetic flux density or the angular 
velocity reach the standard value at n c = 5 x 10 10 cm -3 when 



either B 3s 1.50 fiG or fi c > 2.57 x 10" 



yr 



1 at the stage 



of n c = 5 x 10 4 cm 3 . The latter condition is equivalent to 
|V X v\ > 2.49 x 10 3 kms" 1 pc _1 . 

The standard magnetic flux density is approximately a half 
of the critical one, as mentioned in subsection 4.4. The latter is 
evaluated to be 



B c 



= 2ttVGE 



V 8ivpci 



(19) 
(20) 



at the limit of the geometrically thin self-gravitationally bound gas 
disk iNakano & Nakamuralll97^) . Also, the standard angular ve- 
locity is approximately a half of the critical one. The critical angu- 
lar velocity is defined so that the centrifugal force balances with the 
gravity. Then it is evaluated to be 



4nGp 



(21) 



for a uniform gas sphere. Thus, equation 1 1 71 implies that either 
the magnetic flux density or the angular velocity is regulated to be 
a half of the critical value. 

Equation I17i also predicts anti-correlation between the mag- 
netic flux density and angular velocity of the first core. In other 
words, only one of the magnetic flux density or the angular veloc- 
ity is close to the standard value. Then we can make a new index, 
the ratio of angular velocity to the magnetic flux density, for identi- 
fying whether the magnetic field dominates over rotation during the 
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cloud collapse. If it is larger than the ratio of the standard values 

fist 



= 0.39 VGc 



1.69 x 10 



"' (i 



yr 



(22) 



(23) 



,0.19 kms" 

the centrifugal force dominates over the magnetic force. Otherwise 
the magnetic force dominates over the centrifugal force. This anal- 
ysis suggests that there exist two types of first core: magnetic first 
core and spinning first core. We discuss the difference between 
them in Paper II. 

We shall apply the above discussion to LI 544, the prestellar 
core, of which rotation and magnetic field have been measured. The 
rot ation velocity i s eva luated to be 0.09 km s -1 at r = 15000 AU 



b^JOhash^yd. 
Willi amseTaiT 



< 19991) and 0.14 kms" 1 at r = 7000 AU by 



1999). These velocity gradients correspond to 



th e line-of- s ight magnetic field is evaluated to be +11 ± 2 //G 
bv lCrutcher & Troland fcOOOh . Combining these values, we obtain 
Q/B = 1.1 x 10" 7 yr -1 /xG" 1 and 3.8 x 10" 7 yr -1 /^G -1 . If we 
take account of uncertainty of the observed values, the magnetic 
force dominates over the centrifugal force only marginally. 

It should be noted that lCrutcher et alH 120041) derived a much 
stronger magnetic field (~ 140 /j,G) for LI 544 from linear polar- 
ization of the dust emission. They derived the value under the as- 
sumption that the randomness of the magnetic field can be ascribed 
to turbulent motion. If the magnetic field is as strong as 140 /iG 
at the distance of 10000 AU, the magnetic force should dominate 
over the centrifugal force. However, their method gives a magnetic 
field an order of magnitude stronger compared with the values de- 
rived by the Zeeman effect. Possible systematic errors should be 
examined. 

Next, we discus s the speed of dyn amical collapse in the 
molecular cloud core. Aika wa et a! . (2001) discussed the possibil- 
ity of deriving the collapse speed from the chemical abundance in 
the prestellar core LI 544. They computed chemical evolution in 
a molecular cloud core, assuming that the density evolution is the 
same as that of the Larson-Penston similarity solution, or by a fac- 
tor / slower. They concluded that the observed chemical anomaly 
in LI 544 is consistent with the model based on the Larson-Penston 
similarity solution from comparison with the slow collapse models 
of / = 3 and 10. The model of / = 3 is supposed to mimic a molecu- 
lar cloud core of which collapse is slowed down owing to rotation, 
magnetic field, or turbulence. Our simulation has shown that the 
slowing by magnetic field and rotation is appreciably smaller. The 
slow-down factor is evaluated to be / = 1.22 for model BS and 1.93 
for model CL. The small slow-down factor makes the chemical di- 
agnosis harder. 

Finally, we discuss the effect of the ambipolar diffusion. The 
evolution of magnetically subcritical cloud including the ambipo- 
lar diffusion has been investigated by Basu & Mouschovious (1994, 
1995a, 1995b) under the disk approximation. They showed that the 
magnetically supercritical core is formed in the subcritical cloud 
for ambipolar diffusion after 10-20 freefall time passed. Once the 
supercritical core is formed, the magnetic field is hardly extracted 
from the core, because the ambipolar diffusion is much slower than 
the freefall (Basu & Mouschovias 1994). Thus, our ideal MHD ap- 
proximation is valid since our model cloud is supercritical from the 
initial stage (see Table 2). The ambipolar diffusion may have an im- 
portant role in the initially subcritical cloud and after the formation 
of the dense core (n c > 10 11 cm -3 ). 
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Table 1. Model Parameters 



parameter values 

a 0, IO" 3 , 5xl0" 3 , 0.01, 0.05, 0.1, 0.5, 1, 2, 3 

u) 0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6 

A v 10- 3 , 0.01, 0.1, 0.2, 0.3 

n c ,o 5 X 10 2 cm- 3 , 5 X 10 4 cm- 3 , 5 X 10 6 cm" 3 



Table 2. Parameters and Initial Conditions for Typical Models 



Model a u> A z A v n Cj o -B zc ,o (pG) r2 c ,o (io 8 yr _1 ) M(M Q ) L (io 5 AU) tf(io 6 yr) M/Mb,, 



A 


AS 


0.01 


0.01 


0.1 


0.01 


5 X 10 2 


0.295 


1.26 


12.2 


6.92 


5.96 


13.2 




AL 


0.01 


0.01 


0.1 


0.2 


5 x 10 2 


0.295 


1.26 


12.2 


6.92 


5.99 


13.2 


B 


BS 


0.1 


0.01 


0.1 


0.01 


5 x 10 2 


0.931 


1.26 


12.5 


6.82 


5.92 


4.2 


BL 


0.1 


0.01 


0.1 


0.2 


5 x 10 2 


0.931 


1.26 


12.5 


6.82 


5.95 


4.2 




CS 


0.01 


0.5 


0.1 


0.01 


5 x 10 2 


0.295 


63.1 


20.6 


5.94 


5.35 


14.4 


C 


CL 


0.01 


0.5 


0.1 


0.2 


5 x 10 2 


0.295 


63.1 


20.6 


5.94 


5.38 


14.4 




DS 


1 


0.5 


0.1 


0.01 


5 x 10 2 


2.95 


63.1 


28.7 


5.71 


4.69 


1.4 


D 


DL 


1 


0.5 


0.1 


0.2 


5 x 10 2 


2.95 


63.1 


28.7 


5.71 


4.69 


1.4 
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Figure 1. The density (grey-scale and contour) and velocity distributions (arrows) in model AS [(a, oj, A v )=(0.01, 0.01, 0.01)]. Panels (a) through (f) are 
snapshots at the stages of (a) n c = 5.5 X 10 2 cm- 3 (I = 1,2), (b) 5.3 X 10 3 cm- 3 (I = 2,3), (c) 5.1 X 10 4 c m - 3 (I = 3,4), (d) 9.7 X 10 6 cm" 3 
(I = 6, 7), (e) 8.3 X 10 8 cm" 3 (i = 9, 10), and (f) 7.6 X 10 10 cm" 3 (i = 12, 13), where I denotes the level of subgrid. The level of the subgrid is shown 
in the upper left corner of each solid square which denotes the outer boundary of the subgrid. The upper and lower panels show the cross sections at 2=0 and 
y=0 planes, respectively. The grey-scale is different for each panel as shown above the panels. Thick contours denote the density of n = l/10n c , where n c 
means density at the center. The elapsed time, density at the center and arrow scale are denoted in each panel. 
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Figure 3. The oblateness (upper panels) and axis ratio (lower panels) are plotted against the elapsed time (left panels) and central density (right panels) for 
models AS, AL, BS, BL, CS, CL, DS and DL. The relation of £ a r cx f^J is also plotted in panel (d) for comparison. 
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Figure 4. The square root of the ratio of the magnetic pressure to the gas pressure normalized by the initial value, (-Bzc/Bzc,o)/( n c/wc,o) 1//2 . is 
plotted against the central density, n c , in the left panel. The ratio of the angular velocity to the magnetic flux density normalized by the initial value, 
(f! c /f! Cj o)/ (Bzc/Bz C fi), is plotted against the central density in the right panel. The relation (ric/^co) 1 ^ 6 is a l so plotted in the left panel. 




Figure 5. The magnetic field lines are shown from a bird's eye view for model AS. Panels (a) through (f) denote the same epoch as those of Fig. 0(a) — (f), 
respectively. The density (gray scale) and velocity (arrows) in z = plane are also shown on the bottom. The level of subgrid and central density are shown 
in the upper section of each panel. Each frame denotes a cube with side lengths of (a) 1.7 X 10 5 AU, (b) 1.7 X 10 5 AU, (c) 8.4 X 10 4 AU, (d) 1.1 X 10 4 AU, 
(e) 1.3 X 10 3 AU, and (f) 1.6 X 10 2 AU, respectively. 
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5"'] t= 5 . 97 jhTytl p = 6 Qe+ 1 Ic rr 





50 100 150 



Figure 6. The c 

panels) and y 

(I = 3, 4), (c) 7.9 X 10 6 cm -3 (I = 6, 7), and (d) 6.0 X 10 10 cm -3 (I = 12, 13), The contours, arrows, and notation have the same meaning as in Fig.fTI 



: density (gray scale and contours) and velocity distributions (arrows) for model BL [(a, u>, A v )=(0.l, 0.01, 0.2)] are shown on the z = (upper 
(lower panels) planes. Panels (a) through (d) are snapshots at the stages of (a) n c = 8.2 X 10 3 cm -3 (I = 2, 3), (b) 5.6 X 10 4 cm~ 3 
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Figure 8. The density (grey scale and contours) and velocity distributions (arrows) for model CS [(a, uj, A v )=(0.01, 0.5, 0.01)] are shown on the z = (upper 
panels) and y = (lower panels) planes. Panels (a) through (d) are snapshots at the stages of (a) n c = 5.2 X 10 3 cm -3 (I = 2, 3), (b) 6.5 X 10 4 cm -3 
(I = 3, 4), (c) 5.7 X 10 6 cm -3 (/ = 6, 7), and (d) 8.3 X 10 10 cm -3 (I = 12, 13), The contours, arrows, and notation have the same meaning as in Fig.fTI 




Figure 9. The magnetic field lines at the same epoch of Fig. [8| Each frame denotes a cube with side lengths of (a) 1.7 X 10 5 AU, (b) 8.4 X 10 4 AU, (c) 
1.1 X 10 4 AU, and (d) 1.6 X 10 2 AU, respectively. The gray scale and arrows have the same meaning as in Fig. 151 
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Figure 10. The density (grey scale and contours) and velocity distributions (arrows) for model DL [(a, u>, A v )=(\, 0.5, 0.2)] are shown on the 2 = (upper 
panels) and y = (lower panels) planes. Panels (a) through (d) are snapshots at the stages of (a) n c = 5.6 X 10 3 cm -3 (( = 2, 3), (b) 6.8 X 10 4 cm -3 
(I = 3, 4), (c) 5.1 X 10 6 cm -3 (/ = 6, 7), and (d) 5.3 X 10 10 cm -3 (I = 12, 13), The contours, arrows, and notation have the same meaning as in Fig.fTI 



(a) L=3 nc =5.6x1 3 (cm- 3 ) (b) L=4 ric=6.8xl0 4 (cm" 3 ) (c) L=7 it =5.1x1 s (cm- 3 ) (d)L=l3 ric =5.3x 1 10 (cm" 3 ) 




Figure 11. The magnetic field lines at the same epoch of Fig. fTol Each frame denotes a cube with side lengths of (a) 1.7 X 10 5 AU, (b) 8.4 X 10 4 AU, (c) 
1.1 X 10 4 AU, and (d) 1.6 X 10 2 AU, respectively. The gray scale and arrows have the same meaning as in Fig. 151 
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Figure 12. The evolutions of the magnetic flux density and angular velocity at cloud center. The lower axis indicates the square root of the magnetic pressure 
(B zc / V 87r) normalized by the square root of the thermal pressure (a/ cf p c ) and the left ordinate does the angular speed (f2 c ) normalized by the freefall 
timescale (y/AnGp c ), respectively. The upper axis indicates the parameter, a. The symbols, *, o, o, and a, represent at n c = 5 X 10 2 cm" 3 (initial state), 
5 x 10 4 cm- 3 , 5 x 10 6 cm , and 5 X 10 cm , respectively. Each line denotes the evolutional path from the initial state (n c ,o = 5 X 10 cm i ) to 
the end of the isothermal phase (n c = 5 X 10 10 cm -3 ). The characters 'A', 'B', 'C and 'D' denote the model shown in Table HI The thick curve denotes 
the magnetic flux and spin relation, f! 2 /[(0.2) 2 X AnGpc] + B 2 c /[(0.36) 2 X 87rc 2 p c ] = 1 [see Equation 1171 1. The thin line indicates the relation of 
f2 c / (AirGpc) 1 /" 2 oc Bzc/ (Stc^Pc) 1 / 2 . The models with no magnetic field are shown inside the upper left box, while those with no rotation are shown inside 
the lower right box. 



